Multiplicative noise in the longitudinal mode dynamics of a bulk semiconductor laser 
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We analyze theoretically and experimentally the influence of current noise on the longitudinal 
mode hopping dynamics of a bulk semiconductor laser. It is shown that the mean residence times 
on each mode have different sensitivity to external noise added to the bias current. In particular, 
an increase of the noise level enhances the residence time on the longitudinal mode that dominates 
at low current, evidencing the multiplicative nature of the stochastic process. A two-mode rate 
equation model for semiconductor laser is able to reproduce the experimental findings. Under a 
suitable separation of the involved time scales, the model can be reduced to a ID bistable potential 
system with a multiplicative stochastic term related to the current noise strength. The reduced 
model clarifies the infiuence of the different noise sources on the hopping dynamics. 

PACS numbers: 05.40.-a,42.65.Sf,42.55.Px 



I. INTRODUCTION 

Fluctuations and noise are inherent to the behavior 
of any physical system. Their ubiquity and impact on 
technological applications have stimulated since long ago 
the study of their effects in different branches of science 
(statistical physics, mathematics) and technology (com- 
munications engineering) . 

Noise is usually perceived as a source of degradation of 
the properties of any system, leading to relatively small 
stochastic variations of a given magnitude around its 
deterministic (noise free) value. Such a point of view 
is rooted on the study of systems close to equilibrium, 
whose dynamics can be linearized around the determin- 
istic steady state. However, it does not hold in systems 
far from equilibrium, where nonlinear corrections to the 
dynamics must be taken into account. In these systems, 
noise may not only be the responsible for the observed 
behavior, but it may also make possible new states and 
behaviors that do not exist in the noise- free limit |lj . Ex- 
amples of such behaviors are, for instance, the enhance- 
ment of the decay time of a metastable state (noise en- 
hanced stability) 0, Q , the synchronization with a weak 
periodic input signal (stochastic resonance) or the ap- 
pearance of a periodic output (coherence resonance) [5| . 
In systems with spatial degrees of freedom, noise may 
lead for instance to the formation of convective struc- 
tures (noise-sustained structures) 6]. 

The effects of noise depend on whether it is additive 
or multiplicative. Additive noise is present in all real, 
non isolated systems where the environment acts as a 
thermal bath. Moreover, even in an isolated system all 
interaction processes exhibit some degree of stochastic 
fluctuations that lead to internal noise in the system. 
This effect may be further enhanced by adding noise to 
the control parameters of the system. In general, the 



fluctuation-dissipation theorem implies that noise in a 
nonlinear system may have both an additive and a multi- 
plicative component. The effects of both kinds of noise in 
nonlinear systems have been thoroughly studied from the 
theoretical point of view y| , although experimental stud- 
ies of the effects of multiplicative noise are more scarce. 
In particular, it has been shown that the characteristic 
signatures of either type of noise can be profoundly mod- 
ified by the presence of even a weak component of the 
other kind |3 1^ ■ Moreover it was also found that multi- 
plicative (parametric) noise induces a shift in the critical 
mean value of the parameter that controls an instability 

ilia El. 

In this paper we analyze theoretically and experimen- 
tally the influence of multiplicative noise on the mode 
hopping dynamics of a bistable semiconductor laser. 
Since the pumping current enters into the laser field equa- 
tions in a nonlinear way, adding current noise may lead 
to multiplicative effects. We find that the residence times 
of each mode are strongly affected, and that increasing 
the noise added to the current is equivalent — from the 
point of view of residence times — to a lowering of the 
bias current. A theoretical model for the system under 
study is presented and analyzed, with the aim of un- 
derstanding the observable consequences of the imposed 
fluctuations and the role of laser parameters. A better 
insight is obtained by simplification of this model to a 
bistable potential system. This allows us to discuss in 
detail the effects of current noise on the mode hopping 
dynamics. The results are in good agreement with the 
experimental observations carried out for a bulk edge- 
emitting semiconductor laser. 

The outline of the paper is the following. In Sec. II 
we introduce the rate equations for a two-mode semi- 
conductor laser and we derive from them a reduced ID 
Langevin model that describes the hopping dynamics of 
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the freely operating laser. In Sec. Ill we discuss in de- 
tail the reduced model and, in particular, the effect of 
external fluctuations of the injection current. The limit 
case of a slowly fluctuating current is analyzed as well 
and the results are compared with simulation of the rate 
equations. Sec. IV is devoted to the presentation of the 
experimental results. We draw our conclusions in Sec. V. 



II. THEORETICAL ANALYSIS 

A. Rate equations 

Our starting point is a stochastic rate-equation model 
for a semiconductor laser that may operate in two modes. 
Both modes interact with a single carrier density that 
provides the necessary gain. The two modes have very 
similar linear gains, provided that their wavelengths are 
almost equal and they are close to the gain peak. If 
E± denote the complex modal amplitudes, N the carrier 
density and J(t) the injection current, the model can be 
written as 

1 



E+ =-[(l + za)g+-l]i?+ + v/2^VV^+ (1) 



E- ^-[{l + ia)g^-l]E^ + ^/2D:;N^^ (2) 



N ^ ^[Jit) - N - g+\E+\^ - g_\E_\ 



(3) 



where a is the linewidth enhancement factor |12|. The 
modal gains read 



9± = 



N±e{N- Nc) 
l + s|£:±|^ 



(4) 



where e determines the difference in differential gain 
among the two modes while Nc defines the carrier den- 
sity where the unsaturated modal gains are equal. The 
parameters s and c are respectively the self- and cross- 
saturation coefficients. The ^± are two independent, 
complex white noise processes with zero mean ((C±(t)) = 
0) and unit variance i{^i{t)^*{t')) = S,j5{t - t')) that 
model spontaneous emission. Eqs. IjllSI) have to be in- 
terpreted in Ito sense , thus the average power spon- 
taneously emitted in each mode at any time is given by 
4.DspN. 

The deterministic version of Eqs. H1I3() admits four dif- 
ferent steady state solutions: the trivial solution E± = 0, 
two single-mode solutions — E+ ^ 0, E- ~ and 
viceversa — and a solution where both modes are las- 
ing, E± 7^ 0. The single- mode solution where only E^ 



{E- ) is lasing lacks physical sense for bias currents below 
J+ (J-) given by 



J± 



l±e 



(5) 



The trivial solution E± = is the only stable solution 
for bias currents J < min(J^, J^); for definiteness, we 
consider Nc > 1, hence J_ < 1 < J+. Upon increas- 
ing J, the trivial solution looses stability and the system 
switches to the solution E^ ^ 0, E+ = at the laser 
threshold, J = J-. Further increasing the current, the 
sequence of bifurcations depends quite strongly on the 
parameters s, c and Nc- For c > s, this solution may co- 
exist with the solution E^ ^ 0, E- = or even with the 
solution E^ 7^ 0, E^ ^ (in this case, only for a small 
current range). Finally, the solution E^ ^ 0, E^ = 
prevails. 



B. Reduction to an effective model 

In order to better assess the effects of current noise on 
the modal dynamics, we next reduce the rate-equation 
description to a bistable ID system. In the first place, 
we introduce the amplitude-phase coordinates for each 
mode, 

E± — p±expiilj±. (6) 
Using this standard transformations (see we obtain 



N 



2D,r,N 



9± 



P± + V^D.pN e 



-g± + ^2DspN i 
7[J(i) ~N- g+p^. 



(7) 



9-p'-] 



Since the modal phases do not influence the evolution of 
the modal amplitudes and carrier density, we can disre- 
gard them without loss of generality. It is convenient to 
perform a further change to "cylindrical" coordinates. 



p+ = r cos (j), 



P- 



(8) 



In these new variables, is the total power emitted by 
the laser, and (p determines how this power is partitioned 
among the two modes: ^ = ((/> = 7r/2) corresponds to 
emission in mode -I- (-) only, and intermediate values give 
different power to each mode. 
Using again Ref. we obtain 
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N{1 



5r^ 



'2(j))+e{N- Nc) (l + (Tr2 + 5r'^) cos 20 



(l + crr2) - {Sr^cos2(l)Y 
sin 2<j) NSr'^ cos 2(/) + e (A^ - iVc) (l + or^) 



- 1 



2DspN 



N = 



J-N -r 



{1 + crr'^ f - {Sr^ cos2(j)f r^t(in2(f> 
N [1 + crr^ + 6r'^ cos^ 20) + s {N — Nc) (l + OT^ + (Jr^) cos 20 



(1 



((Sr^ cos 20) 



+ v/2S;PV^, , (9) 
(10) 
(11) 



where we have defined 

The parameter a thus represents the gain saturation in- 
duced by the total power in the laser, while S — given by 
the difference between the coefficients of cross- and self- 
saturation — describes the reduction in gain saturation 
due to partitioning of the power among the two modes. 
We also wish to remark that, although we worked in the 
Ito formalism, this form of the equations would be the 



same also in the Stratonovich representations jl4| . 

The general expression just obtained for the dynam- 
ics is too involved to allow for a detailed analysis of the 
effects of noise. In order to simplify the theoretical anal- 
ysis of this dynamical system, we further assume that: 
(i) The asymmetry of the modal gains is very small, i. 

Nc ^ 1, e <^ 1, d <^ 1] (ii) the laser operates close to 
threshold, so that the saturation term is small, crr^ ^ 1. 
In this limit, we find that, to lowest order in the small 
terms, the dynamics is governed by 



- { N-1- Nar^ 
2 V 

sin 20 



2DspN 



N = ■j{j-N-r^N) , 



(13) 

(14) 
(15) 



so the dynamics of the total intensity and of the carrier 
density N decouple from those of the relative phase 0. 
Thus is driven by the other two variables. 

Moreover, the time-scale for the evolution of the rel- 
ative phase is of second order in the small quantities, 
while that of r and N is of the first order. Since we are 
mainly interested in time scales long enough for the relax- 
ation oscillations of the total power and carrier density 
to be totally damped, and not in the transient dynamics, 
we can consider that r and N have reached the vicin- 
ity of their steady state. By neglecting the spontaneous 
emission term in (|13|) . we have that 



iv.iv„.i±^. (16) 

I + aJ 1 + cr 



A more accurate approximation would be to solve the 
Langevin equations for r and N and to insert the so- 
lution into the equation for in order to describe how 
the fluctuations of the former affect the dynamics of the 
latter. For simplicity, we disregard this issue under the 
assumption that the noise in (|13|l and (|15|) is so weak 



that they simply "follow" J as prescribed by Ijltill . For a 
time-dependent current this approximation is valid only 
if J does not change too fast. For example, in the case of 
an Orstein-Uhlenbeck process that will be considered be- 
low one must require the corresponding correlation time 
T to be longer than the typical relaxation time of the to- 
tal intensity. We shall see that this condition is generally 
met in our experimental setups. 

In this representation the dynamics can be geometri- 
cally visualized as follows. The motion is constrained 
along a manifold (approximated as portion of a circle of 
radius tq) connecting the fixed points. Radial fluctua- 
tions (i.e. fluctuations in the total intensity output) are 
completely neglected. The hopping dynamics is thus ef- 
fectively one-dimensional and described by the "phase" 
variable 0, which determines how the total power is par- 
titioned among the modes. 

Altogether, the above calculation yields the reduced 
model 



a cos 20 -I- 6 



sin 20 



2D^ 
tan 20 



V^^^ (17) 
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where we have defined the parameters 



a = NoSrl = 



1 + cr 



(J-1) 



ea 



b = e{No~N,) = ^—{J-J^) 
1 + cr 



(l + a)(J-l) 



and, for later convenience, we have introduced 



.L = 



(1 + a)N, - 1 



(18) 
(19) 

(20) 
(21) 



Notice that the strength of fluctuations depends on 
the current. At this stage, one may get rid of one of the 
three parameters. For example, one may rescale time t — > 
t/a;£,^ ■v/fl'^^ so that the only independent parameters 
are b/a and D^/a. However, since we wish to understand 
the role of the phenomenological parameters in relation 
with the physical quantities we stick to the form (|17|l . 
This choice is especially useful when current fluctuations 
will be taken into account. 

To conclude this Section, we point out that the 
same equation (fT7|) has been derived by Willemsen et 
al. [l^ [l0| to describe polarization switches in Vertical 
Cavity Surface-Emitting Lasers (VCSELs). The starting 
point of their derivation is the San Miguel-Feng-Moloney 
model ilTj] . We point out that the physical meaning of the 
variable cj) is different from here as it represents the po- 
larization angle of emitted light. Thus the potential min- 
ima correspond to the two orthogonal linearly-polarized 
directions. Also, Nagler et al. |18lll9l | have derived a one- 
dimensional Langevin equation for the mode intensity in 
VCSEL starting from a rate equation model for the field 
intensities in the two polarization directions. Indeed, it 
can be shown that, upon the change of variable sin0 = p 
and redefinition of parameters, rnodel H17|l transforms 
into eq. (19) and (23) of Ref. [ll|. This suggests that, 
upon a suitable reinterpretation of variables and parame- 
ters, many of the results presented henceforth may apply 
also to the dynamics of VCSELs. 



III. ANALYSIS OF THE REDUCED MODEL 

In this Section we analyze some features of the model 
in the case in which J is a constant. In this case, (|17|) 
is a one-dimensional Langevin equation with additive 
noise originating from spontaneous emission. The force 
in Eq. (|17|l can be derived from the potential 



-cos2(?!)- D0lnsin2(/« (22) 



The stationary distribution is thus straightforwardly 
computed as 



P(0) = Po exp(-C/(0,i?^)/i?^) 



(23) 



with Po being the normalization constant. Notice that 
the potential depends explicitly on Z?^ (i.e. on Dgp) 
through the repulsive logarithmic term in (|22|l . This term 
is usually very small for weak noise except at the extrema 
(j) = 0,7r/2 where it diverges logarithmically. As a con- 
sequence, P vanishes linearly there. Physically, this cor- 
responds to the fact that, due to spontaneous emission, 
the modes are never completely switched off. 



A. Bistability and effect of the injection current 

The most important property that can be immediately 
drawn from the form of U is that there exists a range 
of current values J3 < J < J4 for which the system is 
bistable. Within this region, the system has three sta- 
tionary solutions, (j)± (stable) and 00 (unstable). To com- 
pute them, let us neglect the term proportional to Z?^ 
in the definition of U: (f>o ~ (l/2)acos(— 6/a), (p- ~ 0, 
0+ ~ tt/2. The bistability domain is evaluated by the 
condition |6| < a yielding the approximate expressions 



J3 = 



eaJs + S 



J4 - 



ea 



eaJs — S 

ea — d 
-I- 00 



if ea > S, 
otherwise. 



(24) 

Those last expressions have been checked to be a good 
approximation of the exact values computed from bifur- 
cation analysis of the deterministic rate equations in the 
limit of small e and 6. The effect of the term proportional 
to in 122|l is to make the bistability domain dependent 
on the noise intensity, i.e. J^{D^) < J < J4{D^). By 
computing numerically the stationary points of U{4), D^) 
we have found that the width of this interval of current 
values is reduced upon increasing with respect to the 
deterministic case. The current value J controls the sym- 
metry of the potential through the term proportional to 
b. Within the bistable region, at the critical value Js, b 
vanishes and the effective potential is symmetric under 
the transformation — > 7r/2 — 0. In this situation the 
hopping between the two modes occurs at the same rate. 
For weak noise, <ti a, we can easily estimate the two 
potential barriers as 



AU± = 



(a ±6)2 



(25) 



In the neighbourhood of the J ~ Jg, i.e. for {J — Jg) ^ 
(Js - 1): 



S , , (5 ± 2ea , , ^ 

kJs-\) + , A J-Js) . (26) 



[\ + a 



8(1 + a) 



Actually, also the noise strength depends on J: ex- 
panding to second order around J, the definition of Z?^, 
Eq. |(in|), we get 



D4,{Js) 



1 - ^ \ -(J - Js) + 



(1+^)' 



(H-aJ,)2(J,-l): 



(l + aJ,)(J, ^1) 



(27) 
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However, it turns out that this correction to the noise 
is pretty small for the parameters we choose. More pre- 
cisely, one can estimate that, for a given change of J, the 
relative change of MJ± is roughly a factor 10 larger than 
the one of Z?^ 21]. We thus neglect the correction terms 
in lEH) and let = D^{Js). 

In other words, we assume that changing J only affects 
the barriers through formula H26|) . With this simplifica- 
tion the corresponding residence times are given by |20| 



8ttD„ 



TgCxp 



■ exp {^AU±/D4 
5±2ea J - J. 



1- 



(28) 



Here is the residence time at the symmetry point. 
Notice that the prefactor depends on the noise strength 
to leading order. Formula H28(l enlightens the role of 
the asymmetry parameters in determining the residence 
times. Indeed, it predicts that the residence times should 
depend exponentially on the injection current with differ- 
ent rates. In particular, we see that increasing J around 
Js may lead to an increase of both T± if 2ea < d or to an 
increase of accompained by a decrease of r_ if instead 
2ea > 6. To discriminate which of the two cases is of rel- 
evance one must thus know the values of the parameters 
for the laser at hand. This can only be accomplished by 
comparing with experimental data. For the laser at hand 
(see below), the measurements indicate that the case of 
interest is 2ea > S. 



B. Comparison with the simulation of the rate 
equations 

Before proceeding further, we check the accuracy of 
the reduction by comparing with simulations of the rate 
equations. For defineteness, we let e = 0.1, s = 1.0, 
Nc = 1.1, 7 = 0.01 and c = 1.1 and perform differ- 
ent runs for different values of Dgp and for J some 10% 
above threshold. With the above choice, S = 0.05 and 
a — 1.05, meaning that the timescale for the slow vari- 
able (j) is estimated to be about one order of magnitude 
larger than those of the fast ones. We also set a = since, 
as explained above, we expect that the phase dynamics 
should not matter. The largest part of the simulations 
were performed with Euler method with time steps 0.01- 
0.02. For comparison, some checks with Heun method Q 
have also been carried on. Within the statistical accu- 
racy, the results are found to be insensitive the the choice 
of the algorithm. 

We first evaluated the stationary distribution P{(j>) 
as evaluated from the data through the relation (f> — 
atan(|£'+|/|i?_|). In Fig. ^ we compare the expression 
of U with — In P. The simulation data reported there 
correspond to the current value J = 1.197 which we em- 
pirically found to yield an almost symmetric shape of P. 
This value only differs by 2% from the one estimated from 



formula (|21|l yielding Js = 1.19523.... We also compared 
the numerically measured residence times T± with the 
prediction of the reduced model (see again Fig. . Since 
we work with values of the noise that are not extremely 
small, we employed the quadrature formula j2fl] 



1 



P{y) dxdy 



(29) 



This functional form accurately follows the numerical 
data (within then 10% or so). In the inset of Fig. we 
show that the cumulative distribution of residence times 
which is Poissonian as expected from the reduced model. 
Another series of simulations (not reported) for an asym- 
metric case, J = 1.18 yield comparable results. 




FIG. 1: (Color online) Comparison between the reduced 
model and the simulation of the rate equations for J = 1.197, 
Dsp — 7 X 10^^ (the other parameter values are given in 
the text), (a) ~lnP evaluated from simulation compared 
with the analytical result; (b) Average residence times as a 
function of 1/Dsp. The line is obtained from the quadrature 
formula 1291 . Inset: Cumulative distributions of residence 
times for Dsp = 7 x 10"''. The lines are the expected Poisson 
distribution with the same average. 



C. Including current fluctuations 

We include the fluctuation of the injected current by 
letting J — > J + 5J{t). The following considerations 
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holds for an arbitrary time dependencence of 5 J under 
the limitations necessary to eliminate Eq. H13|) . For de- 
fineteness, we may focus on the case in which SJ is a 
Ornstein-Uhlenbeck process with zero average {6J{t)) = 
and correlation time r: 



T 



2D, 



that means 



(SJit)SJ{0)) = L>,/exp(-|i|/T) 



(30) 



(31) 



This choice is suitable to model the finite-bandwith of 
the noise generator employed in the experimental setup. 
Notice that r and the variance of fluctuations Dj = (<5 J^) 
can be fixed independentely. 

The effect of a time-dependent current is twofold. 
First, the potential becomes fluctuating: the a, b and 
coefficients become time-dependent quantities and the 
potential barriers AU±{t) change accordingly. For weak 
noise, and close to the symmetry point SJ <^ {Js — 1) 
they are computed from approximated formula (|26|l with 
J — Js replaced by 5J{t), 



AU±it) 



8(1 -ha) 



{ Js - 1) + 



S±2ea 
8(1 -ha) 



6J{t) . (32) 



The barriers' height depends linearly on 6 J to leading 
order. Obviously, this last expression makes sense only 
when the fluctuating term is sub-threshold i.e. whenever 
the system is bistable (a large enough SJ could always 
occur making the potential single- well) . The second ef- 
fect is on the additive noise strength. Since depends 
on SJ the amplitude of spontaneous emission is renor- 
malized. Indeed, since is S— correlated, the process 
y/2D^{J) ^0 can be replaced by y^2{D^) where the 
average is over the fluctuations of the variable SJ. Using 
the expansion H27|) we get 



{D^) = D^iJs) 



1 



(1 + ^)' 



(l + aJ,)2(J, -1): 



■D, 



(33) 



The coefficient in front of Dj is strictly positive mean- 
ing that current fluctuations always enhance spontaneous 
ones. However, as argued in Section IlII Al for the case of 
static changes of J, it turns out that the relative change 
AU± is larger than that of if Dj ^ {Js — 1)^ (see 
however again note HJ). Therefore, in the following we 
can safely set (D^) ~ D^{Js). 

Although expression (|32|l is sufficient to draw some 
conclusion than can be experimentally tested, it is useful 
to write down also the full Langevin equation associated 
with the problem. Putting al the terms together we find 
that, to first order in 5 J, equation (|17(l transforms to 



4> = -U'{4>,D^)-V'{4>,D^)SJ + .j2D~^i^ (34) 
where 



-V'{c^,D^) 



^,rr- — r f (5 COS 20 + ea) sin 2(^ (35) 
2(1 + a) V ^ J -r \ J 



(for simplicity we neglected again the dependence of 
on the current fluctuations). The multiplicative term can 
be thus derived from the "potential" 



v{4>) = 



ea 



16(1 + a) 



CDS': 



4(1 + a) 



cos2(^. (36) 



For an arbitrary choice of the parameters, V has a differ- 
ent symmetry with respect to U meaning that the effec- 
tive amplitude of multiplicative noise is different within 
the two potential wells. If this difference is large enough, 
current fluctuation will remove the degeneracy between 
the two stationary solutions. 

Non-Markovian equations of the form 1)34(1 have been 
thoroughly studied in the literature (see e.g. [2^12^ 12^ 
(2^ and references therein). Although their full analytical 
solution for arbitrary r is not generally feasible, several 
approximate results can be provided in some limits. In 
the following Section we will discuss the case which is of 
experimental importance. 



D. The kinetic limit 

Altogether, the mode switching can be seen as an acti- 
vated escape over fluctuating barriers given by 1)32(1 . The 
statistical properties of the latter process is controlled 
by the current fluctuations. To assess the nature of the 
stochastic process at hand, it is important to discuss the 
relevant time scales. In particular, one should compare 
the relaxation time within the wells with both r and 
the residence times T±. If Te. < t we are in the colored 
noise case. An estimate of Tr is the inverse of the cur- 
vature at (/)o that is approximatively given by 1/a. For 
example, with the set parameters chosen in Sec. lIII Bl one 
finds Tji ~ 200 in model units. 

In bulk semiconductor lasers the residence times are 
generally much larger than Tji (the switching time be- 
tween the two states). Typically, ~ 5 — 10ns while 
residence times may range between 1 and 100 [is. The 
noise correlation time can be somehow tuned but the 
noise generator limit t to be larger of 100 ns (8.8 Mhz 
is the maximal bandwidth used in the experiment, see 
Section Hvl below) . Moreover, the frequency of the re- 
laxation oscillations of the total power is typically above 
1 GHz, and its damping occurs on a ns time scale, in 
accordance with the assumptions made in the theoretical 
analysis. 

At least to a first approximation, we can thus consider 
the limit of large r. This justifies a further reduction of 
the problem to a kinetic description which amounts to 
neglect intrawell motion and reduce to a rate model de- 
scribing the statistical transitions in terms of transition 
rates. If we consider r as a time scale of the external driv- 
ing we can follow the terminology of Ref. and refer 
to this situation as the "semiadiabatic" limit of Eq. 1(34(1 . 
We will now consider separately two limit cases. 
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1. Slow barrier, rare hops: Tr <C t ^ T± 

This corresponds to the situation in which spontaneous 
emission noise is very weak. In this case the residence 
time is basicahy the shortest escape time, which in turn 
correspond to the lowest value of the barrier (the noise is 
approximatively constant in the current range considered 
henceforth). For the case of interest, S < 2ea we can use 
(|32|l to infer that the minimal values of AU± should be 
attained for 6 J cx ^^/Dj respectively. This yields 



K 



2ea±S 



(37) 



where K is a. suitable numerical constant. The ratio of 
residence times rj is thus exponential in the noise RMS: 



T+ 

7f- = CXp 



- K- 



26 



1 + (J D, 



(38) 



Notice that S controls the asymmetry level: if (5 ^ 2ea 
the two residence times decrease at approximatively the 
same rate. This prediction is verified in the simulations 
reported below and also in the experiment. 



2. Slow barrier, frequent hops: Tr <C T± <C t 

This corresponds to the adiabatic limit in which the 
time scale of the external driving is slower than the in- 
trinsic dynamics of the system ^2ff\. The results of this 
subsection are thus not of direct relevance for interpret- 
ing the experimental results reported here. However, we 
discuss also this regime for completeness and to empha- 
size the differences with respect to the previous case. 

To a first approximation we can here treat current fiuc- 
tuations in a parametric way. The switch time will be the 
average of escape times over the distribution of barrier 
fluctuations, i.e. {T±)sj. Using again the approximation 
(|32|l . since the variable 5 J is Gaussian, using the identity 
(exp/3z) = exp(/32(z2^/2) one gets 



T± ~ T^exp 



2((5±2eCT)2 



Dj 



(39) 



This means that in this limit the residence times increase 
exponentially with the noise variance (i.e. exponentially 
with the square of RMS). In other words, the fluctuations 
always increase the eflFective barriers albeit asymmetri- 
cally. 



E. Comparison with the rate equations in the 
presence of current noise 

The above analysis shows that the effect of multiplica- 
tive noise may strongly depend on the actual parameters 
of the model, in particular of the ratio between e and 
6. Accordingly, the simmetry-breaking effects can vary 



considerably depending on the parameter 6, see Eq. (|38|l . 
To further clarify this dependence and as a check of the 
prediction, we performed some simulations of the rate 
equations (|1I3|I together with 13U|) . We investigated the 
effect of changing the value of the cross-saturation coeffi- 
cient c, keeping the other parameters fixed the same used 
in Section lill Bl The DC value of the current J has been 
again choosen empirically to yield an almost symmetric 
distribution P{(j)) for Dj = 0. 

We set T = 4 X 10"* that correspond roughly to the 
experimental values. In Fig. |21we plot — InP for c = 1.1 
and c = 1.3 corresponding to 6 = 0.05 and 6 = 0.15 
respectively. In the second case the multiplicative effect 
of the noise is stronger and leads to a sizeable distortion of 
the curve while for small asymmetries the curves basically 
mantain their symmetry. 




FIG. 2: (Color online) Simulation of the rate equation with 
current fluctuations: — In P for different amplitudes of im- 
posed noise Dj; for c = 1.1, J = 1.197, P^p = 7 x 10"® (a) 



and c = 1.3, J = 1.194, Dsp = 1.5 x 10' 



(b). 



The same type of behavior can be observed in the de- 
pendence of the residence times on the current fluctu- 
ations. Fig. 13^ illustrates how for c = 1.1 {S — 0.05) 
the residence times decrease at approximatively the same 
rate upon increasing current fluctuations and the ratio 77 
remains very close to 1. On the contrary, upon increas- 
ing the cross-saturation coefficient to c — 1.3 {S = 0.15) 
a sizeable symmetry breaking (Fig. |31d): one of the two 
times remains almost constant while the other decreases. 
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This is in agreement with formulae (|37|l and (|38(l . Indeed, 
the ratio between the decay rates of t] obtained from the 
fits (insets of Fig. |2l is about 40, which is roughly a fac- 
tor 2 of the corresponding value as computed from H38|) 
with the simulation parameters at hand. 




Cun-ent fluctuation (RMS) 




0.02 0.04 

Current fluctuation (RJMS) 



FIG. 3: (Color online) Simulation of the rate equation with 
current fluctuations: Residence times T+ (squares) and T_ 
(circles) for c — 1.1 (a) and c — 1.3 (b). The parameters are 
the same as in Fig. |21 The insets report the ratio of the two 
times and the dashed line is an exponential fit, see expression 



IV. EXPERIMENTAL RESULTS 
A. Experimental set-up 

The experimental setup is similar to the one de- 
scribed in [23, 113 ■ Several edge-emitting semiconduc- 
tor laser were tested: two Hitachi Hip 1400 lasing at 
840 nm and three Sharp LT021 MD lasing at 780 nm. 
These lasers are GaAlAs double-hetero-structure with a 
bulk active region and cleaved, uncoated facets. The 
wavelength separation between consecutive longitudinal- 
modes is around 0.3 nm, and the laser emission occurs 
in a single-transverse mode. The laser package temper- 
ature is stabilized up to 0.01°C and the laser current is 
controlled with a very stable (up to 1 fiA) power sup- 
ply. The laser is optically isolated from the rest of the 



set-up by means of an optical diode that avoid spurious 
back reflection. The total laser emission is detected by 
an Avalanche Photodiode detector (APD) (DC-1.5 GHz 
bandwidth) while the time-averaged optical spectrum is 
measured by an Agilent 86140B spectrum analyzer. Indi- 
vidual longitudinal mode detection is obtained by send- 
ing part of the laser output into a monochromator (res- 
olution 0.5 A) with two independent output slits. The 
monochromator can be set in order to have at the two ex- 
its two different longitudinal modes emitted by the laser. 
Each longitudinal mode intensity is then monitored by 
two APD detectors (DC-1.5 GHz bandwidth) placed be- 
yond the slits. The outputs from these detectors are 
monitored on a LeCroy 7200A (500 MHz analogue band- 
width, 1 GS/s). The power spectra of the signals can be 
monitored using an Agilent E4403B spectrum analyzer. 



B. Longitudinal mode-hopping 

Two control parameters can be directly varied in this 
system: the pumping current J and the temperature of 
the laser substrate Tsub- As described in [53, H^l, the 
optical spectra of the lasers tested show a single longi- 
tudinal mode emission for most of the parameter values 
(J, Tsub)- However, there exist small regions in the pa- 
rameter space where the lasers are bistable and exhibit 
mode-hopping between two longitudinal modes. A basic 
explanation of the destabilization of the leading longi- 
tudinal mode is that, fixing the Tsub and increasing J 
the wavelength of the cavity resonance increases due to 
Joule heating of the semiconductor medium. The gain 
curve peak shifts towards longer wavelength as well, but 
at a larger rate. Eventually, the dominating longitudinal 
mode looses its stability in favour of a longer wavelength 
one, which has become more resonant with the gain peak. 
The same happens when, for fixed J, T^ub is increased. 
This type of transition among lasing modes is well-known 
[T^ and occurs sharply in the parameter space. 

Observation of this regime of stochastic mode-hopping 
has been reported in several papers [27L I29I l30l | . Its most 
relevant features are: (i) the total emitted intensity re- 
mains almost constant though each longitudinal mode is 
switching on and off. In other words, the anti-correlation 
of the two modal intensities is very strong (better than 
—0.95 in our case); (ii) the distribution of residence times 
in a mode, defined as the interval between a switching- 
on and a switching-off event, is Poissonian; (iii) a sweep 
of the pump current across this parameters' region, re- 
veals the existence of an hysteresis cycle for the modal 
emission which is a signature of bistability between the 
two longitudinal modes. It is thus clear that these fea- 
tures are straightforwardly explained by the theoretical 
analysis presented above. 

The study of the fluctuations of the laser emission can 
be carried out by monitoring the modal intensities nor- 
malized to the total intensity output: to/ = ij/itot and 
TO// = iii/itot- The mode labeled with / is the one 
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FIG. 4: Experimental probability distribution of the normal- 
ized modal amplitudes (see text). The histogram is repre- 
sented in logarithmic scale with a grey scale code from white 
(low probability) to black (high probability); J — 95.9mA, 
r,„b = 20.4°C. 



with larger wavelength. These are related to p± defined 
in the theoretical analysis, by mjji = p^/r^. Then, 
it is useful to represent the state of the system in the 
phase space (to/-^/^, to//-^/^). Hence, the variable used 
in the theory is simply the angle that this vector forms 
with the horizontal axis. The strong anticorrelation of 
the modal intensities implies that the modulus of this 
vector is almost constant and equal to one. In Fig. 01 
we plot the the probability distribution function in the 
space (m/^/^, 771//^/^). The maxima of the distribution 
correspond to the two stable modes involved in the hop- 
ping. It is evident that the relevant dynamics occur on 
a portion of the unit circle as assumed in the simplified 
theory. 

The probability distribution functions of the angular 
variable as obtained from the experiment are reported 
in Fig. for different values of the pumping current. For 
every we have included in the histograms all the points 
along the radial direction. Two peaks located around 
(f) = 0.2 and — 1.4 are found, corresponding to the ac- 
tivation of the mode / and // respectively. In Fig.|3^ we 
plot the distribution of for the case of Fig. ^ where the 
modal emission is symmetric, J ~ J^. The measurements 
show that upon changing the pumping current an asym- 
metric situation settles down, where the laser emission 
on one of the two modes is favored. In particular, an de- 
crease (resp. increase) of J with respect the value of Fig. 
0] enhances the stability of the mode labeled by // (resp. 
/), see Fig. ^p^c respectively. In all cases, the logarithm 
of the distribution can be fitted with the function (see 
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FIG. 5: The logarithm of the experimental probability distri- 
bution function of the variable (squares) for three different 
values of the current, Ts„(, = 20.4°C. The line is the fitting 
with the function given in the text, (a) J = 95.9myl, fit- 
ting parameters A = 108.5, B = 0, D = 4.6, C = -10.6; (b) 
J = 95.1mA, fitting parameters A = 85, B = -18.7, D = 
4.9, C = —8.1; (c) J = 96.6mA, fitting parameters A = 
90, B = 13.5, D = 2.5, C = -7.2. 



T//, on the pumping current is shown in Fig. El As pre- 
dicted by formula both times depend exponentially 
on J and the two curves have different absolute values of 
the slopes (remember the remark at the end of Sec. lIII All . 
Incidentally, this implies that the quantity {Tj +Tjj)/2 
increases with the current, meaning that hops become 
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FIG. 6: Average residence times for mode / (triangles) and 
// (squares) as a function of pumping current, Tsut='20.4°C . 
Inset: the ratio Tj/Tji. 



more and more rare. Moreover, the ratio 77 = Tj/Tn 
increases exponentially with the current as shown in the 
inset of Fig. This behaviour is also consistent with 
formula if^ . 



C. Effect of noise addition onto the pumping 
current 

In a semiconductor laser driven by a constant {DC) 
current, the main source of noise is the spontaneous emis- 
sion within the semiconductor medium. This noise source 
cannot be varied directly through the control parameters. 
Fluctuations can be added into the system externally by 
adding electrical noise on the pumping current. This sig- 
nal employed in our setup has zero-mean, a bandwidth 
ranging from 100 Hz to 8.8 MHz and it is AC-coupled to 
the laser bias current. 

In Fig. [7^ we show how the average residence times 
T/ // change upon increasing the current noise level 
around the static value Jg- The times are affected in a 
strongly asymmetric way: while T// is almost unchanged, 
T/ decreases exponentially. In Fig. [3d we plot the prob- 
ability distribution function of the angular variable <f) for 
different values of the pump noise. The experimental ev- 
idence confirms that the asymmetry of the probability 
distribution increases upon increasing the noise level. It 
is interesting to remark that an increase of noise level is 
equivalent, from the point of view of the residence times, 
to a decrease of the DC pumping current. This experi- 
mental result confirms that the imposed noise does affect 
the dynamics in a multiplicative way. The qualitative 
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FIG. 7: Effect of external noise added on the pump current. 
The measurements are performed on a laser sample different 
from the one used for Figs. I4I5I |S| requiring slightely mod- 
ified values for the control parameters to obtain the same 
behaviors, (a) The average residence times Ti (stars) and Tn 
(squares) as a function of the RMS level of the noise added 
to the pumping current, Tsub = 20.5°C, J=92.0 mA. (b) The 
corresponding logarithm of probability distribution function 
of the variable (j) for different levels of current noise (noise 
increases from darl; trace to grey trace). 



agreement between the experimental data and the simu- 
lations reported in Figs.|2jD and|2|3 supports the validity 
of our modelling. The exponential dependence of the 
ratio rj on the RMS of the fluctuation predicted by the 
analysis of the reduced model, Eq. is indeed ob- 

served in the experiment (see Fig. 0. This is another 
indication that confirms that the kinetic description suf- 
fices to capture the mode-hopping dependence on noise 
level. 
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FIG. 8: The ratio of residence times — Ti/Tu as a function 
of the RMS level of the noise added to the pumping current 
(parameters are the same of Fig. |7|l. 

V. CONCLUSIONS 

In this paper we have explored experimentally and 
theoretically the effects of external current noise on 
the mode-hopping dynamics of a bistable semiconductor 
laser. We have shown that the residence times of each 
mode are strongly affected, bearing the typical signatures 
of multiplicative noise. We have investigate a theoretical 
model based on a rate-equation description, where the 
bias current enters parametrically into the evolution of 
the modal amplitudes, hence the multiplicative charac- 



ter of its fluctuations. Numerical simulations of the rate 
equations are in good qualitative agreement with the ex- 
perimental observations. Moreover, the reduction of this 
model to a ID Langevin equation describing activated 
escape over a fluctuating barrier has allowed us to draw 
some predictions (e.g. the dependence of residence times 
on noise strength) and to better understand the role of 
the physical parameters. Depending on their values, the 
fluctuating part of the effective potential may have a dif- 
ferent parity with respect to the static one. This explains 
the observation that imposed fluctuations have different 
effects on the hopping rates (symmetry breaking). 

The reduced model and many of the results presented 
here could be useful to describe also polarization switch- 
ing in VCSELs driven by external noise. Indeed, exper- 
imental data show strong similarities between this phe- 
nomenon 1^ and the longitudinal mode dynamics. 
On the theoretical side, this analogy is supported by the 
fact that the polarization dynamics in VCSELs is de- 
scribed by models that are very similar to the one dis- 
cussed here 0, 0, . The inclusion of current noise 
effects along the same line above reported should be fea- 
sible once the typical time scales and parameters for such 
a class of lasers are evaluated. This work, that we plan 
to undertake in the future, will provide a common theo- 
retical basis to the stochastic dynamics of bistable semi- 
conductor lasers. 
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